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ABSTRACT 

The abundance ratios between key elements such as iron and a-process elements carry a 
wealth of information on the star formation history (SFH) of galaxies. So far, simple chemical 
evolution models have linked [a/Fe] with the SFH timescale, correlating large abundance 
ratios with short-lived SFH. The incorporation of full spectral fitting to the analysis of stellar 
populations allows for a more quantitative constraint between [a/Fe] and the SFH. In this 
letter, we provide, for the first time, an empirical correlation between [a /Fe] (measured from 
spectral indices) and the SFH (determined via a non-parametric spectral-fitting method). We 
offer an empirical version of the iconic outline of iThomas et al.l (l2005l [20101) . relating star 
formation timescale with galaxy mass, although our results suggest, in contrast, a significant 
population of old ( ^ 10 Gyr) stars even for the lowest mass ellipticals (Ms 3 x IO^'^Mq). 
In addition, the abundance ratio is found to be strongly correlated with the time to build up 
the stellar component, showing that the highest [a/Fe] ( ^ +0.2) are attained by galaxies 
with the shortest half-mass formation time ( ^ 2 Gyr), or equivalently, with the smallest ( <, 
40%) fraction of populations younger than 10 Gyr These observational results support the 
standard hypothesis that star formation incorporates the Fe-enriched interstellar medium into 
stars, lowering the high abundance ratio of the old populations. 
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1 INTRODUCTION 

The use of spectral indices has allowed a relatively successful 
description of the unresolved stellar populations of galaxies in 
■ terms of three luminosity- weighted parameters: the age, the global 
' metallicity, and the [a/Fe] abundance ratio. The advent of new 
tools of analysis (e.g. spectral fitting), the increasing refinement 
of the population synthesis m odels and libraries (e.g. MILES, 
ISanchez-Blazquez et alj l2006al) and the access to vast spectro- 



scopic surveys (e.g. SDSS, lYorketalj|2000l) allows for a fresh 



look at well established standards. One of the current paradigms 
involves the existence of a relationship between [a/Fe] and the 
star formation timescale. On the basis of this relationship lies the 
assumption that the stellar abundance ratios quantify the contri- 
bution of different stellar chronometers (e.g. supernovae) to the 
chemical enrichment of the subsequent generations of stars. Core- 
collapse supernovae, triggered soon after the onset of star forma- 
tion, release mostly Mg and other a-process elements, whereas 
type la supernovae, with a more complex time delay, contribute 
large amounts of Fe (see e.g. lOreggio & Renzini|[T98 3l:lMaoz et al.l 



1201 iHMannucci et al.ll2005h . Hence the [a /Fe] abundance ratio is 
a useful tool to tr ack star formati on timescales At ^ IGyr (see e.g. 
Matteucci & Tomambe 1987 : Thomas et alJI 19991 : iFerreras & SilM 

boosh . 

By using simple chemical evolution models, iThomas et al.l 

( 12005 !)" established a theoretical relationship between the mass of 
a galaxy and its SFH, described by a Gaussian function. This 
Gaussian curve is characterized by the mean, given by the SSP- 
equivalent age, and the FWHM, related to the [a /Fe] abundance 
ratios, such that [a/Fe] ~ 0.2 - 0.17 log(FWHM/Gyr). There- 
fore a high abundance ratio would result from a concentrated burst 
of star formation, while solar, or sub-solar ratios would reveal a 
more e xtended SFH. The app lication of this rule to observational 
data by IThomas et al. I 1 I2OO5I) and Thomas et al. (2010) produced 
an iconic outline (figure 9 in Thomas et al. 2010) in which mas- 
sive galaxies show older and narrower stellar age distributions com- 
pared to their less massive counterparts. 

There are tools currently available to produce more detailed 
SFHs than the Gaussian curves described above. The use of spec- 



©2011 RAS 



2 de la Rosa et al. 



tral fitting (e.g. Icid Fernandes et ^l2005h allows us to probe the 
full range of populations, giving a more robust estimate of mass- 
weighted properties. Furthermore, our methodology follows a non- 
parametric approach, with the goal of establishing a distribution of 
ages that better maps the spectroscopic data. These estimates are 
less prone to the biases inherent in the use of luminosity-weighted 
properties (when using SSPs) or in the parametric assumptions 
(e.g. when using exponentially decaying SFHs). The analysis of 
the abundance ratios follows instead the traditional method based 
on line indices. We emphasize that both properties (SFHs and 
[a /Fe] ratios) are therefore extracted independently. In the future, 
the development of a-enhanced populatio n synthesis models (e.g. 
ICervantes et"ai]|2007l; ICoelho et alj|2007h will open up the possi- 
bility of a fully-consistent analysis of the SFHs and the abundance 
ratios from a spectral fitting point of viewQ 

The goal of the present study is to revise the link between 
the SFH and \a / Fe] in elliptical galaxies, im proving over the in- 
sightful outline of lThomas et alj ( l2005l , |20 1^ by using a combina- 
tion of spectral fitting, to determine non-parametric SFHs, and line 
strength analysis to derive abundance ratios. 



2 SAMPLE AND ANALYSIS 

Our methodology involves the analysis of the optical spectra to 
constrain the properties of the underlying stellar populations. The 
sample is extracted fr om the ^40,000-strong SPIDER catalogue 
of early-type galaxies ("L a Barbera et alj |2010') (hereafter Paper I), 
which covers a redshift range of 0.05 to 0.095 and absolute mag- 
nitudes brighter than Mr = -20. The spectroscopic data come from 
the seve nth data rele ase (DR7) of the Sloan Digital Sky Survey 
(SDSS) (i Abazaiian e t al. 2009). We take a subsample with a S/N 
high enough for a detailed analysis of the stellar populations. We 
select a set of 9,999 ETGs, belonging to the upper quartile of the 
distribution in S/N, as measured in the ^-band (per-pixel S/Ng > 
18.3). 

Stellar populations are characterized by the following infor- 
mation: (I) mass-weighted average ages ((age) m) and metallicities, 
(II) mass-weighted star formation histories (SFHs) and (III) aver- 
age [a/Fe] abundance ratios. Several techniques and models are 
available for the extraction of this information. In an ongoing study 
(de la Rosa et al. in preparation) the performance of several tech- 
niques and models has been evaluated for SDSS spectra of similar 
quality. The methodology chosen for the present study, as described 
below, results from that performance comparison. 

I) Average ages and metallicities are obtai ned through spectral 
fitting using the STARLIGHT synthesis code J cid Fernandes et akl 
12005) to find the optimal mixture of single stellar populations 
(SSPs). For the present stud y, the basis SSPs co rrespond to 76 
solar-scaled MILES models ( Vazdekis et al. l l2O10h. These models 
have a Kroupa Universal Initial Mass Function ( lKroupj|200ll) and 
span a range of 4 metallicities, from Z/Zq = 1/50 to 1.6, and 19 
different ages, from 0.5 to 12.6 Gyr. These models are b ased on 
the MILEI^ stellar library dSanchez-Blazguez et al.ll2006al) . which 
combines both a dense coverage of the stellar atmospheric pa- 
rameters and a relatively high and nearly constant spectral res- 



Table 1. Uncertainties on Parameters 



Parameter 




typical rms 


SFH age-bars 


percent 


13.4 


(age>M 


Gyrs 


2.7 


<[Z/H])m 


dex 


0.10 


[a/Fe] 


dex 


0.10 


log(Mdyn/M0) 


dex 


0.07 



olution (FWHM=2.3A). The fitting interval spans from 3600 to 
7350A, with emission lines and bad pixels being masked out. 
The extinction due to foreground dust is modelled with the CCM- 
law jCardelli. Clavton & Mathijl 19891) . The standard output of the 
spectral fitting consists of the fractional contribution of each SSP to 
the flux at a normalization wavelength Ao=4020A, converted into 
a stellar mass fraction by the M/L of the given SSP at Ao, as pro- 
vided by the models. The detailed population mixture of a galaxy 
can thus be condensed into a single parameter, the mass-weighted 
average age, (age) m- The mass-weighted metallicity, ([Z/H])^/, is 
defined in a similar way. 

II) The use of spectral fitting with solar-scaled MILES models 
provides the fractional contribution of each SSP to the total flux, 
allowing for the reconstruction of the SFH. Note that STARLIGHT 
does not assume any parametric form, describing the SFH by an 
arbitrary mixture of the basis SSPs. In order to account for the range 
of redshifts of the sample, we add the corresponding lookback time 
to the population ages, ranging from 0.67 Gyr (z=0.05) to 1.24 Gyr 
(z=0.095). Throughout this paper we assume a standard ACDM 
cosmology with Q,m ~ 0.3 and Ho=70 km/s/Mpc. 



m) The 
iThomas et al 



t/Fe 



^ As a robustness c heck, we have derived galaxy SFHs also with the 
ICervantes et al ] i2007l) a-enhanced SSP models, finding good agreement 
with the results based on the widely used, solar scaled, MILES models. 

^ www.iac.es/proyecto/miles 



is measured with the recent models of 
1201 II) . which provide absorption-line index values 
for a variety of ages and element abundance ratios. Our procedure 
fixes the age to its SSP-equivalent value, as obtained from spec- 
tral indices H^ and [MgFe]'. Then, the Mgb and Fe3 = (Fe4383 -l- 
Fe5270 + Fe5335)/3 indices are compared to the a-enhanced model 
predictions to obtain the (SSP-equivalent) [cc/Fe] . To this effect, 
all measured indices are preliminarly corrected for both velocity 
and instrumental dispersion broadening. The correction functions 
are estimated by convolving a variety of SSP-model spectra with 
broadening functions of different widths. 

Dynamical masses are computed with the formula Mdyn,n ~ 
[K(n)Recr^]/ G, where n i s the S ersic index and the function K(n) 
is taken from lBertin et ai] ( I2OO2I) . This assumes that mass follows 
light in a galaxy, and accounts for structural non-homology of 
galaxy light profiles, as parametrized by the Sersic n. The three 
relevant parameters, n. Re and a, are taken from Paper I. 



2.1 Error Estimation 

To estimate uncertainties on stellar population parameters, we have 
resorted to the 2,283 ETGs in the SPIDER sample with repeated 
spectral observations from the SDSS database (see Paper I). Af- 
ter selecting for pairs of repeated observations with S/Ng >18.3 
(i.e. that of the sample used here), a sub-sample of 393 pairs of 
spectra is explored with the same methodology using differences 
in the extracted parameters as an estimate of the expected errors. 
These uncertainties are summarized in table[T] It is worth mention- 
ing that the reported errors refer to single measurements of stel- 
lar population parameters, i.e. they reflect the statistical uncertainty 
inherent in the extraction of SFHs from the noisy spectrum of a 
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Figure 1. Variation of the star formation liistory witli respect to dynamical mass (each panel representing a mass bin as given in table|2)- Sketchy (theoretical) 
SFHs, plotted as grey Gaussian curves, are supe rimposed over the de tailed (empirical) SFH histograms (black filled) obtained through spectral fitting. The 
Gaussian curves are obtained from equation 4 of iThomas et alj f2005h . using our SSP-equivalent (age)^ (+ lookback time) and [ct/Fe] values, with curves 
normalized to span the same area as the histograms. Arrows mark the position of the mass-weighted (age)jv/, obtained from the SFH histograms. Each panel 
also displays Mz<i, i.e. the stellar mass percentage formed after z=l. The grey insets at the upper-left comers are the metallicity histograms, ([Z/H])jv/ , for 
each mass bin, with a dotted line at the solar metallicity. 

Table 2. Stellar population properties as a funcion of M^yn . For each mass bin, en'ors on a given quantity are the root-mean-square of deviations from the 
median value in the bin. 



Bin (logMdyn/M©) 


12.0 


11.6 


11.0 


10.5 


Log Mass interval 


11.8-12.2 


11.5-11.7 


10.9-11.1 


10.3 - 10.7 


Bin members 


1024 


1466 


1408 


690 


(age) A/ (Gyr) 


10.82±1.68 


10.35±1.91 


9.29±2.22 


7.71±2.45 


SSP-equiv. age (Gyr) 


9.07±2.08 


8.74±2.15 


7.97±2.47 


6.61±2.66 


SSP-equiv. [c»/Fe] 


0.17±0.09 


0.15±0.09 


0.11±0.09 


0.07±0.10 



given galaxy. On the other hand, binned quantities in (Figs. 1-2) 
have much smaller uncertainties, resulting from average-stacking 
a large sub-sample of galaxies in each bin. The comparison of re- 
peated measurements also shows no systematics in stellar popula- 
tion parameters as a function of S/N, as shown, for cro, in Paper I. 



3 RESULTS 

3.1 The dependence of SFH with Mdyn 

Mean SFHs are constructed for four mass bins selec ted to match 
the mass intervals of figure 9 in iThomas et all ilOld ). Our results 
for these four mass bins are presented in Table [2] where one can 
see that the age, metallicity, and abundance ratio parameters ac- 
tually increase with stellar mass, consistent with previous studies. 
Figure[T]shows the average distribution of stellar ages for the ETGs 
within each mass bin. We emphasize that constraints of stellar pop- 
ulations from spectroscopic data are more robust when dealing with 
relative age differences. By constructing a stacked age histogram 
for all galaxies within each mass bin, we give a rough estimate of 



the range of stellar ages to be found in galaxies of the correspond- 
ing stellar mass, i.e. an averaged star formation history. Note we use 
seven time intervals only (instead of the original 19 time steps from 
the basis SSPs) in order to achieve more robust results (see dis- 
cussion in lCid Femandes et al.ll200^ . The histograms combine all 
possible metallicities within a given age bin. It is worth mentioning 
that the absence of the youngest populations (t ^ 2 Gyr), evident in 
the depleted left side of the histograms, reflects the limit imposed 
by the lookback time related to the redshift range of the sample. We 
do not consider this a problem since most of the stellar populations 
have ages Si 5 Gyr. Our extracted SFHs ( filled black his t ogram 
in fig ure [D are compared with the ansatz of IThomas et al.l j2005l , 
bold) (grey Gaussian curve), which relate the mean and variance 
of a Gaussian SFH with the age and [ct /Fe] of the underlying stel- 
lar populations. The mean of the Gaussian SFH is set to our SSP- 
equivalent age with an added lookback time, while the FWHM is 
extracted fro m the [a/Fe] 0.2 - 0.17 log(FWHM/Gyr) rela- 
tion given bv lThomas et al.l j2005l) . 

We have also constructed stacked metallicity histograms for 
all galaxies within each mass bin, displaying them as grey insets in 
figure [T] Histograms show high percentages (60-80 %) of super- 
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Figure 2. Trends between the SFH and [a /Fe] . The top panel shows the 
correlation with the mass fraction in old stars ( 10 Gyr). The insets illus- 
trate individual cases for bins 1 , 5 and 9 of the [a /Fe] distribution (grey 
pentagons). The shaded histogram bins in the insets mark the ^ 10 Gyr 
populations used for the definition of the OLD parameter. The bottom panel 
shows the trend with half-mass time (i.e. the time interval over which 1/2 of 
the total stellar mass is formed). Individual cases of the stellar mass growth 
for bins 1, 5 and 9 are shown in the insets, with the black lines guiding the 
eye for the value at 0.5 of total mass. 



solar metallicities, with a growing contribution (from 6 to 20 %) 
of lowest metallicities ([Z/H]<-0.5) from the highest to the lowest 
galaxy mass bin. This low metallicity trend might be real, or just 
reflect the age-metallicity degeneracy and the fact t hat low-[Z/H] 
stars in the MILES library have super-solar [a /Fe] ('Milone et al.' 
l201lh . To test if this can bias our results, we have re-computed the 
age histograms by excluding galaxies with a significant (> 20%) 
low [Z/H](<-0.5) population. We found no significant variation 
wrt figure [T] (e.g. less than 1% variation in the SFH for age > 
10 Gyr). 

A c omparison between our SFH extraction and the Gaussian 
ansatz of IXhomas et"ai] ( |2005|,|201(]|) . shows an overall agreement, 
although our non-parametric approach gives a more detailed view 



of the underlying SFHs. We emphasize that it is the relative dif- 
ference across mass bins what should be compared and not the 
absolute age values, whose estimate can be significantly affected 
by the adopted methodology (e.g. use of spectral-fitting vs. in- 
dices) and stellar population libraries. Notice also that the cen- 
ter of the Gaussians does not coincide with that of the SFH his- 
tograms, as expected by the fact that luminosity-weighted ages are 
systematica lly younger tha n their mass-w eighted counterparts (e.g. 
ISerra & Tra ger 2007; Rog ers et al.ll2010l) . Both the SFH and the 
Gaussians show the characteristic trend towards younger ages and 
more extended age distributions in lower mass galaxies. 

However, figure [T] goes beyond a simple description of the 
SFH as a Gaussian. The non-parametric approach of STARLIGHT 
allows us to study in more detail the relative change of SFHs 
with respect to galaxy mass. One notices that over the mass range 
probed, all galaxies feature a significant population of old stars. If 
we consider the stellar mass formed at redshifts z> 1 (i.e. older than 
~ 8 Gyr in a local sample) we still get about 40% of stellar mass 
formed even at Mdyn ~ 3 x IO^^'Mq. Furthemore, the most mas- 
sive galaxies (Mdyn ^ lO^^M©) undergo an intense, short-lived 
period of star formation, in which over 50% of the stellar mass is 
formed within a very short period (At ;S 2 Gyr) at redshifts z ^ 4, 
followed by a tail of low-level star formation. After z~ 1 , only 18% 
of the stars are newly formed, perhaps the contribution from minor 
mergers that could explain the size evolution of massive galaxies 
(see e.g. lNaab et al.ll2()09l:lTVuiillo et alj|201 ih . 



3.2 The dependence of SFH with [a /Fe] 

In order to confirm a link between the SFH shape with [a/Fe], 
we bin the sample with respect to the abundance ratios into 9 bins, 
each one containing 1,111 objects. An average SFH is constructed 
for each bin and two parameters are defined to describe the dis- 
tribution of stellar ages. The OLD parameter measures the stellar 
mass fraction contributed by the oldest populations, namely above 
10 Gyr. The upper panel of Fig. |2]shows the increase with [a/Fe] 
of the OLD parameter. A simple linear regression gives for the mass 
fraction in old stars: 

OLD = L397[Q/Fe] + 0.344. (1) 

For clarity, three insets show the average SFHs for bins 1, 5 and 9 
(corresponding to lowest, intermediate, and highest [a/Fe]). The 
distribution of the extracted values of the OLD parameter for galax- 
ies within the same bin has a standard deviation around 25%. The 
OLD parameter in these insets corresponds to the added height of 
the two shaded oldest bins of the histograms. According to this 
result, the drop in the [a/Fe] abundance ratio is motivated by 
a growing fraction of younger populations on top of the preva- 
lent old com ponent. This result supports the standard interpreta- 
tion (see e.g. lSanchez-Blazguez et al]|2006bl) , whereby old stellar 
populations, formed from primeval ingredients, show high [a /Fe] 
abundance ratios, further reduced by subsequent star formation 
episodes, which incorporate the Fe-enriched inter-stellar material 
into stars. 

A complementary view of the same result uses the cumulative 
integration of the SFH to calculate the "half-mass time" (Tjvf/2)> i-e. 
the time needed to form 1/2 of the final stellar mass of the galaxy. 
We note that the initial time is set to a lookback- time of 13 Gyr for 
all galaxies. However, this is an arbitrary choice depending on the 
basis of SSP models used for the spectral fitting. The lower panel 
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of Fig. |2] shows a decreasing [a /Fe] abundance ratio related to a 
growing half-mass time revealing the connection between an 

extended period of star formation and the decrease of the [a /Fe] 
abundance ratio. Similarly to the top panel, the three insets show 
the average cumulative distribution of the SFH for bins 1, 5 and 9. 
A simple regression gives: 



TM/2(Gyr) = -15.3[a/Fe] + 5.2, (2) 

with a scatter around 2.5 Gyr (rms) for galaxies belonging to the 
same bin in [a /Fe] . The fiducial value for solar abundance ratios 
is represented, therefore, by a fraction of 34% in stars older than 
10 Gyr and a half-mass formation time around 5.2 Gyr. 

In general, absolute values of stellar population parameters 
depend significantly on methods and models used to derive them. 
Hence, we warn the reader that equations (1) and (2) apply exactly 
only to the set of models used in the present work, as well as to the 
lookback time scheme adopted here to bring all SFHs into a com- 
mon time-frame. In other words, usin g equations I and 2 all ows one 
to convert the [ct /Fe] measured with lThomas et al. models 
into half-mass time Tm/2 and mass fraction in old ( ^ 10 Gyr) 
stars, as they would be measured with solar-scaled MILES models 
and lookback times added up. 



4 CONCLUSIONS 

This letter presents an independent comparison between the [a /Fe] 
abundance ratios and the star formation history of a volume-limited 
sample of ~ 10* giant early-type galaxies (Mr < —20), ex- 
tracted from the Sloan Digital Sky Survey. Figure [T] shows the 
characteristic trend towards more extended SFHs in lower mass 
ETGs, in agreement with previous work (e.g. iThomas et aLl l2005h 
but with the additional information regarding the distribution of 
stellar ages. We find that all ETGs in our sample have a signifi- 
cant fraction of very old stars (ages ^10 Gyr), at variance with 
the si mplified ansatz of a Gaussian SFH of iThomas et alj ( l2005l 
l20IOh . The fraction of stars fanned after z ^ 1 increases from 
18% at Mdyn ~ IO^^Mq up to 57% for Mdy„ ~ 3 x 10i°Mo. 
This points to the young populations, formed out of a Fe-enriched 
medium, as responsible for the decrease of the initially high 
[g/Fe] abundance ratios of th e old stellar components (see e.g. 
ISanchez-Blazquez et alj2006bh . 

In addition to the SFHs, we used individual absorption lines 
to deter mine the [ct/Fe] abun dance ratios, following the prescrip- 
tions of IThomas et al.l ( l201lh . Figure |2] shows the correlation be- 
tween these two independently estimated properties. We find a 
strong correlation between [a/Fe] and either the mass fraction in 
old ( ^ 10 Gyr) stars (OLD parameter; top panel), or the time to 
form one-half of the current stellar mass of the galaxy (Tm/2, bot- 
tom). Both cases give a clear picture towards enhanced [a/Fe] for 
the SFHs that are more short-lived. Those for which Tm/2 ~ 2 Gyr 
feature highest, super-solar, abundance ratios, with [a/Fe] ^ 0.2. 
Alternatively, an [a/Fe] ^ 0.2 is reached by those galaxies where 
the old populations ( ^ 10 Gyr) constitute more than 60% of the 
stellar mass content. These quantitative results are powerful cali- 
brators of models of galaxy formation. 
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